r.C T '".A I. F^J?C HT r^CTlOII 

T wADUATl SCBOCft 
vr \ 4 CAXIFG&KIA 9Mtt 



NPS-67Bi77111 

NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




AERODYNAMIC STABILIZATION OF GASEOUS DISCHARGES 



0. Biblarz 



Final Report for Period 1 Oct 1976 - 30 Aug 1977 



Approved for public release; distribution unlimited 



FEDDOCS 

D 208.14/2:NPS-67Bi771 11 



NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



Rear Admiral I. W. Linder 
Superintendent 



Jack. R. Borsting 
Provost 



This report reflects work performed during AY 1976-77 while the author 
was on temporary duty at the Department of Aeronautics at the Technion in 
Haifa, Israel. 

The author would like to thank his colleagues and students at the 
Technion for helping make his stay a memorable experience. 



Reproduction of all or part of this report is authorized. 



This report was prepared by: 



UNCLASSIFIED 

SECURITY CLASSIFICATION OF This PAGE (W*on Dm* Kntmrob) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1 4tro*T HUMIC* 2. QOVT ACCCtSIOM NO. 

NPS67B177111 


2. RECIPIENT'S CATALOG NUMBER 


4 Tl TL E find Sub nil*) 

AERODYNAMIC STABILIZATION OF GASEOUS DISCHARGES 


». Tyre or Mtrorr t remoo covimeo 
1 Oct 1976 - 30 Aug 1977 


s. PERFORMING ORG. REPORT NUMBER 


7 AuTHORfi^ 

0, Biblarz 


S CONTRACT or grant numb t*(m) 


S PERFORMING ORGANI 2 A T ) ON NAME AMO AOORESS 

Naval Postgraduate School 
Monterey, California 93940 


10. PROGRAM ELEMENT. PROJECT, TASK 

AREA i WORK UNIT NUMBERS 


n controlling orricc nauc and Aooress 

Naval Postgraduate School 
Monterey, California 93940 


12. REPORT DATE 

1 Nnvpmhp^ 1977 


it. Nuwacr or r age* 

28 


1 4 MOniTORinG AGInCY namC A AOORESVl/ d///«r«ni from Controlling Olllcoj 


it. security class. r«> u»<« r^panj 

UNCLASSIFIED 


ii«. oecl assi fi cation/ downgrading 
SCHEDULE 


1*. DISTRIBUTION STATEMENT (ol t*lm Mepo ft) 

Approved for public release; distribution unlimited 


17 DISTRIBUTION STATEMENT (ot (ho ebotrect on toted In Blook 20, 1/ different from 1 Report) 


»B. SUPPLEMENTARY notes 


19 KEY WORDS (Continue on tor or •• ml So If no coeoery end Identify be bio ok number) 

Glow Discharges 
Stability 
Aerodynamics 
Lasers 



20. ABSTRACT (Continue on r«r«r«« «/*« i/ Atc«««A^ Identity by bio ok mm m b • * ) 

This report covers some aspects of the formulation of the stabilization 
problem by aerodynamic means. Results of experimental work are not included but 
are referred to in appropriate places. The classical solution for ambipolar 
diffusion in tne positive column of the glow discharge must be extended to re- 
flect the influence of current and pressure. The influence of current is given 
by the inclusion of appropriate energy equations and the pressure influence by 
the addition of two-body recombination. It was found that, from dimensional 
analysis, the characteristic lengths for ambipolar and heat diffusion for 



DD , 1473 

(Page 1) 



coition or i nov «i it oatoccTc 

S/N 0102*014* 4*01 I 



UNCLASSIFIED 



i 



ICCliaiT’*' CL AMI P' CAT ION OF TMI* RAC1F :r 



UNCLASSIFIED 

iuCu witv Classification of This p »GE r'*'s#n n*r« Er * t*fi 



helium do not match. A re-formulation of the problem is thus necessary. The 
recombination term does not change the eigenvalue nature of the diffussion 
equation except at the stationary or local equilibrium limit. 



DD Form 1473 
1 Jan 73 

S/N 0102-014-6601 



UNCLASSIFIED 

ii SECURITY CLAtSiriCATlON OF TW»S FAGt/^w Dmf 



TABLE OF CONTENTS 



Section Page 

I. INTRODUCTION 1 

II. SPECIFICS 2 

III. DISCUSSION 3 

IV. CONCLUSIONS 11 

APPENDIX A - PROBLEM REFORMULATION 12 

APPENDIX B - INFLUENCE OF TWO-BODY RECOMBINATION OF SCHOTTKY 
SOLUTION 17 

REFERENCES 21 

LIST OF FIGURES 23 

DISTRIBUTION LIST 24 



I. INTRODUCTION 



In the high energy laser field, the electric discharge convection laser 
(EDCL) offers good efficiency as well as simplicity of handling relative to 
gasdynamic and chemical lasers. Progress towards higher laser outputs from 
EDCLs, however, has been impeded by the collapse of the glow discharge at the 
increased currents and pressures desired. This glow collapse (usually re- 
ferred to as arcing or constriction) represents the practical upper limit of 

the electrical power that may be pumped into the necessary non-equilibrium 

1-4 

processes. It has been shown that aerodynamic stabilization can be used 
very successfully in overcoming the glow-collapse limits. It should be 
pointed out here that the laser application is but one application of the 
relatively new interdisciplinary field of cold plasma chemistry. 

Aerodynamic stabilization is based on two separate effects usually re- 
ferred to as convection (laminar) and mixing. The purely convective effect 
is said to decrease the residence time of the particles through the discharge 
and thereby to counteract unstable mechanism that lead to the glow collapse^. 
The mixing effect, which may be due to vortex flows or to turbulence from 
boundary layers or grids, has been observed to have powerful stabilizing 
characteristics. Convection scales up with flow velocity whereas mixing 
depends on the intensity, spectral distribution, and location of the eddies 
or vortices . 
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II. SPECIFICS 



At the Technion, the author had the opportunity to participate in 

experimental work dealing with vortex-flow stabilization^’^ and with super- 

6 8 

sonic flow stabilization ’ . These techniques complement the work that has 
been underway at the Postgraduate School on turbulence stabilization of 
discharges. Some results of the Technion work have already appeared (see Refs. 
6-11) and more will be forthcoming. 

A significant portion of the author r s time was dedicated to a study of 
theoretical aspects of discharge stabilization. The author was privileged to 
colaborate with Dr. Y. Khait in establishing a perspective regarding some of 
the author's previous results with turbulence stabilization. This effort 
resulted in the writing of a paper titled "Influence of Turbulence on Diffuse 
Electrical Gas Discharges Under Moderate Pressures," which is being submitted 
for publication. Other theoretical work involved a close look at the 
Schottky solution for the positive column in a glow discharge with the aim 
at establishing its use for a stability formulation. This effort was fruit- 
ful on two intermediate aspects of the problem and these are included here 
as Appendixes A and B. 

A proposal was written to the United S tates-Israel Binational Science 
Foundation^ which extensively outlines the progress of the entire effort. 

If awarded, the cooperation of the Technion and the Postgraduate School 
would be formalized (however, the budget would provide neither salaries for 
principal investigators nor equipment for the USA experiments) . 
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III. DISCUSSION 



The classical solution to the ambipolar diffusion of charges in the 
12 13 

positive column * of a glow discharge offers no hints as to how striations 
and constriction, among other important observations, manifest themselves in 
the discharge. As pressure and/or current are increased, instabilities grow 
in the discharge and it is important to have a sufficiently complete model 
of the physics which reflects these effects. In other words, before flow 
effects such as (laminar) convection and turbulence are introduced it is 
necessary to re-examine some assumptions inherent in the Schottky solution. 

A partial list of phenomena not reflected in the Schottky solution may 
be given as follows: 

1) effect of sheaths at the walls 

2) effect of electrodes and intermediate regions 

3) effect of recombination in the gas (pressure and composition) 

4) effects of the energy balance, including radiation and the elevation 
of the electron temperature 

The list can be rather long. We closely examined the energy balance because 
the instabilities are believed to be thermal disturbances (i.e., the tempera- 
ture becomes non-homo geneous) and recombination in the gas because it closely 
reflects the effect of pressure and composition. For simplicity, a two- 
dimensional (infinite parallel plate) geometry was considered, see Fig. 1. 
When E is constant, 7*J = 0 implies that J = J(y) or that 

Li X 

n = n(y) because = 0 ; moreover, V x E = 0 implies that = E(y) . 



The quasineutral equations for steady state may be written as 



10,14 
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Figure 1. Geometry of Positive Column 
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" 3/3y 2T F (nT e ) . = V ‘ an2 


(1) 


- 3/3y [K 3T/3y] = j-E = eny E? 

e L 


(2) 


2 — V f k(T - T) = ep E 2 
DL. eN 2 e e L 


(3) 


p = NkT 


(4) 



with boundary conditions 

at y = 0 3n/3y = 3T/3y = 0 

at y = + a n *= 0 and T = T 
— o 

In the above formulation x is the axial coordinate and y the interplate 
coordinate, n the charge concentration in the ambipolar region, the 

electron temperature and T the gas/ion temperature, N is the concentra- 
tion of the neutrals, the coefficient of ambipolar diffusion, K the 

thermal conductivity, m^ the electron mobility, and a the ionization 

and recombination coefficients respectively, E the axial electric field, 

Li 

v the collision frequency between electrons and neutrals, e the elec- 
tronic charge, p the pressure, and m^ and m^ the mass of the electrons 
and neutral respectively. There are several assumptions inherent in this 
formulation, namely that N>>n , T^>T , and that p = constant • In addition, 

we are using a symmetry condition in the boundary conditions and we are 

12 13 

neglecting of wall space charges. Glow discharge data 9 , in addition, 

indicate that as well as appear to be constant in the discharge at 

least before constricting. The electron energy equation, Eqn. 3, should have to 
reflect the energy dependance of the collision cross sections and, in addition, 
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non-elastic collisions — since we are interested in laser operation. 

However, we shall not deal with Eqns . (3) and (4). 

Using simple kinetic theory to describe the gas temperature dependance 

of the coefficients (D , u , K , etc.) and defining a stationary state 

a e 

where 

v.n* - an*^ = 0 (5) 

i 

or n* = v j/ a 



we can proceed to nondimens ionalize Eqns. (1) and (2). We shall, however, 
not use the plate spacing to render y dimensionless, but introduce two 
characteristic lengths, and , The gas temperature is made non- 

dimensional by the use of a wall temperature T q . The resulting equations 
are: 



, , , /rrl/2 dn N 2 

- d/dy (T — ) = n - n 


(6) 


IT D" 
t - / e ao 

^ " V 2v. T 
V lO o 


(7) 


- d/dy (T 3/2 ~) = nT 


(8) 



/ T K 
1 o o 

I = V e n* u E 2 

eo L 



(9) 



It seems 
glow, the two 



reasonable to expect that in a conventional (unconstricted) 
characteristic lengths be about equal to each other and that 



L^/a a 0(1) and L^/a ~ 0(1) 



( 10 ) 
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This expectation is also required by the eigenvalue nature of Eqns. (6) and 
(8). It may be shown (see Appendix B) , that except for the equilibrium limit, 
Eqn. (6) is an eigenvalue equation, (the classical solution being a cosine 
when no recombination is present and when T « constant) . The nature of 
Eqn. (6) is discussed further at the end of this section. 

In the light of Eqn. (10), it is interesting to look at some actual 
values of and . Particularly extensive information for helium is 

reported in Ref. 16, and this is used in estimating the values shown in 
Fig. 2, which are limited to the range of applicability of the equations 
in Ref. 16. An appropriate form for the coefficients is obtained by intro- 
ducing the gas density N q , 



L 

n 



N 

o 



T (D /N ) 
e a o 

T 

O 1 o 



( 11 ) 



L_ N = (E/N ) 
1 o o 



-1 



T K a 

o o 



e(y N )(y./N ) 
eo o x o 



( 12 ) 



The above forms show "similarity" since LN is a function of E/N only. 

Examining Fig. 2, ic is evident that L^N is not equal to L^N , certainly 

not for large E/N , wheareas L^N follows closely data reported for helium 

12 13 

in the literature (the fact that we are looking at parallel plates and 

the data is for tubes is immaterial). Since in a discharge tube n and T 
share a common diameter, it is necessary to reconcile the discrepancy of the 
two curves . 

14 

It is possible, that the assumption which brings T^ into Eqn. (1) 

and, therefore, into Eqn (11) is invalid. But, even if we were to eliminate 
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HELIUM 



ro 

i 
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LOG (LN) 



the factor T^/T q the equation, the curves would diverge at the higher 

values of E/N . Constriction of the positive column is another possibility 
but in this case would be expected to be greater than for any given 

E/N and just the opposite seems to be true. Uncertainties in the estimate 
of a , the recombination coefficient, would affect the value of but 

this estimate is the one that follows nicely the experimental values. 

The method of reconciliation of the two lengths proposed here is based 
on the alternative that the equations previously written are incomplete and 
that, in particular, in Eqns . (1) and (7) should be replaced by a 

parameter (C^) which makes up for the difference in the two curves . This 
matter is discussed in Appendix A of this report. 

Consider now Eqn. (1) with constant coefficients. We proceed here to 
nondimensionalize it in a conventional way 

- D n" = v. n - ctn^ (13) 

a l 

n T (o) = 0 and n(+ a) = 0 



Let n(o) = n and e = an /v. 

o o 1 

Hence we get 

"2 

- n = n - en 

n T (o) = 0 and n(o) = 1 

n(+ a/ D /v.) = 0 
a l 



(14) 



Now the Shottky solution is one for which e = 0 , namely, 



and 



n = n cos 
o 



(f y/a) 



(15) 

(16) 
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It can be seen therefore that the characteristic length L should be 



L = 




a 

ir/2 



(17) 



On the other hand, for the equilibrium limit e = 1 and the equation 



n 



A 

n - n 



(18) 



is of the "Boundary Layer" type, i.e., it is not an eigenvalue problem, see 
Ref. 14 pp 149-51 for the equivalent problem with three-body recombination. 
The solution to Eqn. 18 in an implicit form is 



y = log 



(n + i ) 1/2 






- 1.317 



V 2 “ ( n + 2 ) 



1 . 1/2 



(19) 



or explicitly 



a 



3 

*T n 
2 o 



1 _ e - (y + 1.317) 

l 7 e - (y + 1-317) 



n 

o 



2 



( 20 ) 



which is indeed of the boundary-layer type. The intermediate cases, 

0 < e < 1 , are represented by elliptic functions, and this matter is dis- 
cussed in Appendix B of this report. 
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IV. CONCLUSIONS 



Some deficiencies of the Schottky solution have been re-examined with 

the intent of establishing a more useful solution for the description of 

the transition to the glow collapse. At least two major, recent attempts in 

the literature (Refs, 16 and 17) base the steady solution for the charge 

profile on the Schottky result and these attempts have not succeeded in pre- 

10 . 

senting a complete formulation of the problem. The Technion approach is 
to base the perturbation on a physically complete solution of the ambipolar 
charge profile. This profile must include the ingredients of constriction 
and of the elevation of the electron temperature. 

As shown in Appendix A, the characteristic length scales for ambipolar 
diffusion and heat diffusion (conduction only) may be reconciled if in 

Eqn. (7) is replaced by C ? where and C ^ is a function of 

E/N (i.e. , the dependance given by the difference of the two curves in Fig. 
2). What appeared to be a bonus, namely, a insight into the mathematical 
description of striations, needs further study in order for it to make more 
physical sense. The eigenvalue nature of the diffusion equation is not 
charged by the recombination term except at the limit when e = 1.0. 
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APPENDIX A - PROBLEM REFORMULATION 

Equation (1) in the text presupposes steady flow along with n = n(y) . 
Let us relax these suppositions and go back to the species flux equations 
which may be written as 

, _ D 

# + V • [n(- u E - V n) ] = n (Al) 

o t e n 



I 7 + v- [n(u. E - -i V n)] = h 
a t In 



(A2) 



Combining and eliminating the electric field we get the ambipolar diffusion 
equation in the usual manner which for constant coefficients is 



or 



3n/3t - D 

a 

3n/3t - D 

a 




2 

v. n - an 
1 



(A3) 

(A4) 



The above assumes two-body recombination as well. Now, using Eqn. (A3) to 
2 

eliminate V n from Eqn. (A2) yields 



8 n 

3t 



+ 



V* nE 

1 - D./D 
1 a 



= n 



(A5) 



But since V*E = 0 in the ambipolar region 



3 n *i 

at 1 - D./D 
i a 



E • V n = n 



(A6) 



or 3n/3t + 



U . 
i 



, - dn , _ on V 2 

1 - D./D I E Lh +E ,a’ # i n ' 



(A7) 



i a/ 
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Let us work with Eqns. (A4) and (A7) in linear form, i.e., neglecting 
recombination. If 



n(x,y,t) = f (y) g(x,t) 
Then, substituting into Eqn. (A4) 


(A8) 


f g t ' °a (f g xx + gf,,) = V i fg " * * * 


(A9) 


d 2 f 3 2 

where f" = — j g = — | and g = 3g/3t 

dy 301 3x 

Dividing by fg both sides of (A9) 

1/5 (8 t - D a - r £ " ‘ \ 


(A10) 


At a given E/N , is constant so we may presume that 

1/g (g - D g ) = constant (or - C ) 

L <3. XX J. 


(All) 


D 

and - — f = constant (or C^) 


(A12) 


with - C, + C_ = v. 

1 2 l 


(A13) 


The constants and C 7 are positive and defined so that C ^ 


V . 
1 



as required by Eqn. (A12) , which now replaces Eqn. (6) since 



Equation (A14) is the new Schottky equation. The constant appears in 

the companion partial differential equation 



g - D g + C g = 0 
6 t a 6 xx 1 & 



(A15) 



Now look at Eq. (A7) . At the centerline 3n/3y = 0 and f(o) = f 



so that 



at y = 0 



3n ^ 111 V- 3 n _ 

3t + \1 - D ± /D a y E L 3 x V i n “ 



(A16) 



Introducing n = fg we obtain 



S- + 



U i 



E t g„ f = v. f g 



’to \1 - D./D / L s x o i o l 
l a / 



(A17) 



We can cancel f in the linear equation above and obtain an equation 
in which no function of y appears 



8 t + \l - bjvj E L 8 x ■ V 



(A 18 ) 



At this point it is worthwhile to eliminate g^ between Eqns. (A18) and 
(A15) and obtain 



°a 8 xx “ °1 8 V i g 1 - D./D / E L g > 

\ i a/ 



P i E L 

°a 8 xx + (1 - D./D) s x ~ v "i ' "l ; 
1 3 . 



g = (V. + C )g 



or 









8„ = 7T 8 



3 xx \ D - D. / °x D 
\ a 1 / a 



(A19) 
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Equation (A19) can be solved with 



g = T(t) e 



ax 



(A20) 



where a is found from 



- 5 ± \/« + 4C./D u E 

a Ji— , 2—S. l L 



D-D. 
a i 



(A21) 



Note that D_^ < D^ and that £ is anyway a positive parameter. 
Equation (A21) indicates that the roots are real and there is no periodic 
behavior in x — that striations are not to be found with this description. 

In order for n to oscillate in x , we must have a situation in which 
new constants C^* and C ^ are introduced such that 



and 



C l* - C 2 * - V. 



g - D g = C * g 
& t a & xx 1 & 



- D f" = - C * f 
a 2 



(A22) 

(A23) 

(A24) 



Here, Eqn. (A24) will result in (non-periodic) hyperbolic functions in y and 
Eqn. (A24) in oscillatory functions in x provided that 



4 r * 
4L 2 

D 

a 




(A25) 



As before we would like to have 

C l* ? V i (C 2* = C l* “ V i ) 

so as to match L to L„ . 

n l 



(A26) 
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The solution to Eqn. (A24) does not represent the usual diffusion 
controlled situation because in order to meet the symmetry condition at the 
centerline we must give up the physically sensible boundary conditions at the 
wall. Perhaps variable coefficients and the inclusion of recombination can 
remedy this problem with the boundary conditions. Or perhaps we should be 
willing to give up axial symmetry and consider some sort of contraction 
asymmetry which would likely have to include some reconbination because of 
the higher pressures. This, however shall not be pursued in this report. 
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APPENDIX B - INFLUENCE OF TWO-BODY RECOMBINATION OF SCHOTTKY SOLUTION 



The diffusion equation given in the main text (Eqn. 14) is not suitable 
for elliptic funciton representation, therefore, a different form is developed 
below. 

The diffusion equation with constant coefficients is (without considera- 
tion of whether C^ shall be used) 

n d^n 2 /'R-n 

- D — 7 T = v. n - an (Bl) 

a , 2 i 
dy 

2 

Defining a stationary state vm* - an* = 0 , or 

n* = \>Ja (b2) 

We can nondimensionalize the equation with 

n = n/n* and y = y/L 



Obtaining 




- 2 
n 



(B3) 



n' (o) = 0 and n=B, 0 <_ S <_ 1 

n(a/L ) = 0 
n 



It is clear that the boundary condition at the centerline make up for the 
disappearance of t in Eqn. (B3) . Now Eqn. (B3) can be integrated once to 
yield 



dn 

dy 





(B4) 
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or 





dn 

yi3 t~2 ~t 

n - y n - p 



3 2 

+ f 6 




(B5) 



The left-hand-side of the above equation can be represented in terms of ellip- 
tic functions. Let 



-3 3-2 2 3 2 

(n - f a - + f fT) = (n - a) (n - b) (n - c) 



therefore 




b = g (0 < g < 1) 



c = 




where a > b > n > c 



(B6) 



(B7) 



With information on a, b, and c we can enter any table for elliptic 
functions (see for example Ref. 18 or 19). Using Byrd and Friedman notation, 
formula 234.00, we have 



where 




A 

dn 

\/(a - h) (b - n) (n - c) 



gF (v?,k) 



g = 



2 



V a - 



ic * 




(B8) 
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and 



<P = sin 



. -1 




(a - c) (b - n) 
(b - c) (a - n) 



F(<p,k) is tabulated. Note that for 3 = 1, we get a double root (a = b) 
and the equation is not an elliptic integral (though its solution can still be 
found from the Tables) . When 3=1 we get the boundary-layer behavior 
described by Eqn. (19) or (20) in the main text. 

Except for 3 = 1 , Eqn. B3 represents an eigenvalue problem. In order 
to check this, we solve for the value of the integral at the wall where 
n = 0 . Thus 



As 3 increases, we obtain progressively larger values until at 3=1 we 
get A = 00 ) . Figure B1 shows this behavior for the eigenvalue. The results 
shown in Fig. B1 proved identical to results from a computer program at the 
Technion (for the combined diffusion and heat dissipation). 

Increasingly higher degrees of recombination flatten out the density 
profile at the center and increase its slope at the walls. The radial elec- 
tric field is thus enhanced. The variable coefficient solution is not antici- 
pated to change the above results in any major way. 




For very small values of 3 , we obtain the cosine solution, namely 



A = (tt/2)^ where L = a/A^“ 

n 
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Figure Bl. Eigenvalue for Increasing Degrees of Recombination 
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